Skin friction in zero-pressure-gradient boundary layers. 
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A global approach leading to a self-consistent solution to the Navier-Stokes-Prandtl equations for 
zero-pressure-gradient boundary layers is presented. It is shown that as Res — s- oo, the dynamically 
defined boundary layer thickness S(x) oc xl In 2 Re x and the skin friction A = ^S>. oc 1/ In 2 5(x). Here 

t w and Uo are the wall shear stress and free stream velocity, respectively. The theory is formulated 
as an expansion in powers of a small dimensionless parameter i^l q in the limit x — >• oo. 



The law of variation of skin friction with Reynolds num- 
ber in turbulent wall flows is one of the oldest riddles 
of physics of turbulence. In addition to difficulties as- 
sociated with a general problem of strong isotropic tur- 
bulence, the presence of solid walls is responsible for ap- 
pearance of two different characteristic velocities. The so 
called friction velocity, reflecting properties of the near- 
wall sublayer, is equal to m 2 = v\ -7^ \waii and the ex- 
ternally prescribed free-stream velocity is denoted as Uo- 
For dimensionless distance to the wall y+ = ^f- = 0(1), 
the dimensionless velocity profile U + (y + ) — U{y + )/u* oc 
is independent upon the Reynolds number. It is the 
interplay of the two characteristic velocities that makes 
derivation of mean velocity distributions in wall flows so 
difficult. 

The first semi-empirical theory of zero-pressure-gradient 
boundary layer (BL) was developed by Prandtl (for de- 
tailed description see Ref.[l]) and references therein) 
who, denning self-similar variables U(x, y) — U{jj^pj), in- 
troduced the boundary layer thickness S(x) depending 
upon distance to the origin x. In addition, from the 
force balance, Prandtl obtained the expression for skin 
friction A oc 4^ oc (rf-) 2 ■ Then, as a stroke of genius, 



based on experimental data on skin friction in pipes and 

channels X p i pe oc Re D 4 available in those early days, he, 
identifying S(x) with the pipe diameter D, proposed the 
differential relation: 



Then substituting into this expression the assumed veloc- 
ity profile ^jf^- oc <f>(sk\) with either algebraic or loga- 
rithmic expressions for the function <j)(z) led to different 
relations for skin frictions reasonably well representing 
experimental data in a restricted range of Reynolds num- 
ber variation. Due to the lack of mathematically sound 
derivation, even today, the question of " correct" velocity 
distribution in fully developed pipes/channel flows is a 
topic of lively scientific discussions. These early theories, 
based on bold assumptions combined with deep physical 
insights, led to helpful expressions widely used in engi- 
neering. 

The modern analysis of pipe/channel and boundary lay- 
ers is typically based on an assumed scaling relation for 
velocity represented in the "inner" and "outer" regions 
of the flow as U(y) = u st f{y + ) and U(y) = U C L — Uog(jj), 
respectively 2 . The parameters u* and uq are correspond- 
ing characteristic velocities. Then, different matching 
conditions applied to the "overlap" region lead to dif- 
ferent shapes of velocity profile U(y). According to ex- 
perimental data 2 , in the limit Ren — > 00, the second 
characteristic velocity uq — > u*. In a recent paper 3 
assuming the logarithmic velocity profile across a zero- 
pressure-gradient boundary layer, developed an asymp- 
totic expansion, leading to the so called Coles-Fernholtz 
relation 4 : 



dS(x) 
dx 



oc( 



SU ' 



giving S(x) oc R ^ Ll , and A oc Re x - the expressions of- 
ten used in modern engineering literature and textbooks. 
A different, more detailed, formulation was based on a a 
force balance following the equations of motion: 
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In Re S ( x ) 
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widely accepted as an accurate large Reynolds number ( 
Reu » 10 5 ) asymptotics. While this work is based on a 
solid mathematical analysis, its starting point, logarith- 
mic profile, is an assumption not following the Navier- 
Stokes equations. 

The theory developed below is based on the following ob- 
servation. It is well known that in the small-parameter- 
lacking homogeneous and isotropic turbulence evaluation 
of moments of velocity fluctuations is an extremely dif- 
ficult and unsolved problem. In this regard the prob- 
lem of statistical properties of wall turbulence is at least 
as difficult. However, in wall flows as Res — > 00, the 
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d5(x) 
dx 







global dimensionless parameters u* /U — > and 
are small and can be used for construction of the well- 
behaved perturbation expansion leading to prediction of 
global properties of wall flows. This is the main goal of 
this paper. 

We consider a flat plate < x < oo and y = 0. The 
freestream velocity of incoming flow is Uo = Uq'i and we 
are to analyze the Navicr-Stokcs -Prandtl equations in 
the boundary layer approximation: 



dU dV 
dx dy 



= 



U— + V— - —{v— +T 

dx dy dy dy x ' v 



(2) 



(3) 



with r XtV = —u x u y and the mean is taken over ensemble 
of fluctuating velocity field u. Interested in the global, 
integral, features of a flow we neglect contribution com- 
ing from extremely thin in the limit Res — > oo viscous 
sublayer < y < y S L — > and as x — > oo we, assum- 
ing self-similarity of velocity profile write: U (x, y) = 
U{x,jfc) = U(x,T]), V = V(x,-gfa) = V{x,rj) and 
T x , y = T Xi y(x,-J^j) = T Xi y(x,r]) where the defined be- 
low width of the boundary layer 5(x) must be found from 
equations of motion. The incompressibility constraint (2) 
gives: 



V(x,y) = - f 
Jo 



v d U(x,y') 
dx 



dy' 



v dU(x,rf) 



dS(x) r> ,dU(x,r]') 



Integrating (3) over the interval < y < oo, and in- 
troducing the 'displacement thickness' 8 we, using (4), 
express the skin friction in terms of the boundary layer 
thickness S: 



<m__ r 

dx 7 



G(x, V )( 



S dU{x,rf) 



Uo dx 



| dS{x)U{x,rj) 



dx Uo 



frr x/2 



(5) 



where ul(x) = v^Ml\o and G(x, V ) = 1 — ^). In 
this work we are interested in the global, integral, fea- 
tures of wall turbulence, like U with the negligible con- 
tribution from the viscous sublayer y/S < y s h/5 —¥ in 
the limit Res — > oo. Therefore, in what follows we seek 
solution in the form 



U(x,y) = u*(x)<p(ri) 
so that the relation (5) becomes: 



(6) 



dx Uo Jo 
d6(x)u*(x) f°° 
J o Jo 



dx Un 



= 172 = *i + «2 = A /2 



(7) 



It will become clear below that in the limit Re x = 
xUo/v — > oo the first term in the right side of (5) is 
negligibly small and therefore A = 2(|^) 2 oc^§^0. 
Based on as yet unknown function S(x), we define an 
averaged- over- the -boundary -layer property *: = 
| Jq ^(x,y)dy = Jq ^(x,r])dr]. Since at the edge of a 
boundary layer y = S(x), the velocity and kinetic energy 
are U = U(x,S(x)) and K = K(x,5(x)), respectively, 
the familiar integral balance equations must be some- 
what modified. For example, integrating the differential 
energy balance equation: 



TJ dK jr dK du c d . dK 

> - -1 — = - -£ + -^{v— + Q) 



dx ' ' dy xy dy 
in the interval < y < S(x) and recalling that 



f s dK [ s dV 

Iv =nwi iV 

we, using an incompressibility constraint, derive: 



f -K(x, 
Jo dx 



y)U(x,y)dy + V(x,6)K(S) 
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Jo 



Txy-^jdy - SE + Q(5) 



(8) 



where Q(5(x)) = w(5)v%(5) = 0{u\). In the limit 6 — > 
oo, V(x, 5) — > and the relation (8) tends to a familiar 
energy balance (see for example Ref.[5]). With 



V(x, 6(x)) = -5^ + £(U(x, S) U{x)) = 

~Tx^ + Tx U ^ < 9 > 

directly following from (4), the relation (8) takes a very 
simple form: 



^(SKU) K{x, S)-^(SU) = jf 4 rj-^dy - 5E + Q(5) 

(10) 

The left side of (10) is the energy balance in the cross- 
section x. As Re x — > oo we can define K (x) so that the 
mean: 



3 



KU=- s jy K(x, y)U(x, y)dy = K 



(x)U(x) (11) 



In the bulk of the flow y 3> y S L dominating 
the integral (11), where the turbulent kinetic energy 
K (x, y) is a weak function of y, the parameter (K(x) — 
K (x))/K(x) 1. 

Expansion. Now, recalling that = u*(x) depends on 
the distance to the origin x, we introduce a definition of 
the boundary layer thickness: 



U(5(x)) - U(x) ocU - U(x) = u^Mtt) 

Uq 

K (x)=u 2 J K ((^) 2 ) 
K(5(x))=u 2 A((^) 2 ) (12) 

As Res — > oo the dimensionless parameter £(x) = 
— > and taking into account that in this limit all 
functional derivatives 



d n U+ 



d 



n Kg 



d£ n ' 

are finite (for a discussion see below), we derive: 

oo 

U(S(x)) - U(x) oc Uq - « u* V a„(J)" 

OO 

K (x) = uij2b n (^r 

n—0 

oo 

^(x))« U ^ Cn( ^ )2n (13) 
where K n (x) is defined as: 



71=0 



It will become clear below that as x — > oo, ^uj^r^ = 
0{K/Uq). Substituting the anzatz (12)-(13) into the 
energy balance (10) and taking into account the estimates 
(15), we, equating the terms of equal powers of small 
parameter u^/Uq obtain in the zeroth order: 

f s BJJ 
J o T xy —dy = 5£ + Q(6)=0(ul) 

Locally, the dissipation rate £ satisfies a well - known 

4 4 

universal scaling relation £(y) = -^-£+{^jr) = -jr£+(y+) 
which has recently been tested experimentally in pipe 
flows 6 - 7 . Thus: 



•/ 

Jo 



5E = ui / £+(y+)d 



r R, 

z / 

Jo 



£+{y+)dy + = 0{ul) 



justifying the above zero-order estimate, provided in the 
limit i?* — >• oo the function £ + decreases with y + rapidly 
enough so that the convergent integral is independent of 
the upper limit. The production term is: 
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dn = 0(ul) 



(r + = r xy /u 2 ) where in the limit Res — > oo the integral 
is dominated by the outer scaling U{y) = Uq — u*g(y/5) 
(Smits (1999)). The flux contribution Q(S) = 0(ul). 
According to this estimate, contributions to the right and 
left side of (10) are balanced separately. This fact be- 
comes clear if we notice that each term in the left side 
of (10) involves x-derivatives which are proportional to 
^| oc (ff^) 2 — > while each contribution to the right side 
is 0{u%) = 0(1). Moreover, since in the interval y > 5, 
all turbulence characteristics sharply decrease to zero, all 
integrals over the y- coordinate in the right side are con- 
vergent so that the outcome is independent neither of the 
upper limit 5 nor The remaining terms are: 



KU = K (x)U(x) = ul(U -u* V ) 4 )x V b n &) 

1=1, U ° rTT^O U ° 

, . (W) 

The relations (14) is a definition of expansion coefficients 
b n . Let us demonstrate that the anzatz (12)-(13), com- 
bined with the energy balance (10), leads the well-known 
empirical relation A oc ln 2 Re ■ Neglecting for a time be- 
ing the contribution i\ to (7) as small (it will be justified 
below) we obtain simple estimates: 



In 



—ul = -4) < -f ; 

dx U§ Uq 

(15) 



dx 



n=0 

oo 
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n=0 
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We see that the first -order in u^/Uq terms are balanced 
only if Co = bo and in the next order we, neglecting the 
asymptotically small 0((f^) 3 ) contributions, obtain the 
differential equation: 



TT ?dul » dS 

Uod-j— = ctiu^ — 
dx dx 



(17) 
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where ct\ — 2— — — < is a constant which must be 
obtained from a full local theory. 

It is easy to see that the expression: 



M* . 2 _ k dS 
U~o> ~h^~5 ; dx~ 



A = 2 (^) 2 = 



(18) 



with k = 8/ct 2 is a solution to (17). Indeed, integrating 
(17) and dividing the outcome by Uq , we obtain: 



A 
2 



AV2 



dx = K^ A 

dx 4V2 



This result shows that the anzatz (12)- (13) with A oc 

is a self-consistent solution to the Navier-Stokes -Prandtl 

equations of motion. 

Setting for a time being all proportionality coefficients 
equal to unity, we introducing So — ^r, Reg = and 

R e x = ^fr solve the the differential equations (17) with 
the result: ilea [(In Res) 2 - 21n-^-] = Re x leading to 



S(x) oc 



In 2 Re x 



(19) 



in the limit Res — > oo. 

The derivation of asymptotic expressions (18)-(19) was 
based on the relation oc (ff^) 2 oc A which is obtained 
from (5) by neglecting the first contribution to the right 
side. With the scaling (6) this is equivalent to omitting i\ 
in (7) as small. Now, based on the derived solution (18)- 
(19) we justify this approximation. Indeed, comparing 
the integrands in the right side of (5) gives 



dx Uo 



dV\ 



In Re x dx 



G 



G 



In 4 Re x 



(20) 



and 



»2 



dS 

dx Uq 



G 



In 3 Re,, 



(21) 



so that as Re T — » oo the ratio 



0. Thus the term i\ 



in (7) is asymptotically negligible. 
Summary and discussion. 

1. In this work the expression for skin friction in the 
zero-pressure-gradient boundary layer has ben derived 
directly from the Navier-Stokes equations. Similar re- 
sults were recently presented in a resent paper 3 . In this 
work, to fit experimental data, the authors used a par- 
ticular forth-order Pade approximant yielding logarith- 
mic velocity profile in a leading order. This immediately 
gave U + (S) rts \nS + and the Coles- Fernholtz relation (1) 



for the skin friction. Since it has not been shown that 
the Pade approximant is an asymptotic solution to the 
Navier-Stokes equations, this theory is essentially semi- 
empirical. 

2. The theory presented here is different. No exper- 
imental information has been used to derive (18)-(19). 
Moreover, while we can say that these expressions are 
consistent with U+ oc lny + , the mathematical proof that 
this is so does not exist. 

3. In arriving in the final result (18)- (19) we assumed 
that all scaling functions </>*(£) are analytic. It is impor- 
tant to understand that in a logarithmic theory one can 
expect singularities when dealing with derivatives over 
spatial coordinates x or y. It is not what has been done 
deriving the expansions (12)-(13). It is clear that the 
derivatives over the functions £(x) = w*(x)/J7o cannot 
be singular by definition: it is impossible to obtain a 
large response of, say, U+ upon infinitesimal variation 
of the friction velocity w*. This simply contradicts all 
conservation laws. 

4. The expression (18)-(19) for friction factor (skin 
friction) was found as a self- consistent solution to the 
Navier-Stokes- Prandtl equations. The full mathematical 
theory leading, in addition to the scaling relations derived 
here, to the amplitudes, remains an unsolved problem. 

5. These results are accurate up to the O(u*/Uo)- 
corrections. 
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